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Abstract 

As a follow-up to a recent study in the spin-0 case [J. Bunao and E. 
A. Galapon, Ann. Phys. 353, 83-106 (2015)], we construct a one-particle 
Time of Arrival (TOA) operator conjugate to a Hamiltonian describing 
a free relativistic spin-1/2 particle in one spatial dimension. Upon trans¬ 
formation in a representation where the Hamiltonian is diagonal, it turns 
out that the constructed operator consists of an operator term T whose 
action is the same as in the spin-0 case, and another operator term To 
which commutes with the Hamiltonian but breaks invariance under par¬ 
ity inversion. If we must impose this symmetry on our TOA operator, 
then we can throw away To so that the TOA operator is just T. 


1 Introduction 

The marriage of relativity and quantum mechanics is a difficult one. Even for 
the case of special relativity, the interpretations can be somewhat troublesome. 
In the regime where the effects of quantum mechanics and special relativity 
are no longer negligible, the particle number in a system can, and will, fluc¬ 
tuate. Roughly speaking, this is because mass is just another form of energy, 
according to special relativity, and energy can fluctuate, according to quantum 
mechanics. Thus, there can be massive particles being created and annihilated 
in a relativistic quantum system. This poses a problem on the interpretation of 
the supposedly one-particle wavefunctions satisfying relativistic wave equations 
(i.e. the Klein-Gordon and Dirac Equations) and the quantum operators rep¬ 
resenting the physical observables of the system. To illustrate, consider a time 
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of arrival (TOA) experiment. Suppose we start with a particle in some initial 
state and let it propagate towards a detector. When a particle arrives at the 
detector, it records the time of arrival. However, how can we be sure that the 
particle that arrived is the same one we started with? This is not an issue in 
the usual quantum mechanics but it is a question that must be addressed in 
the relativistic quantum regime. We can also repeat the same experiment with 
particles starting with the same initial states in order to get a TO A probability 
distribution. How can we then interpret this distribution not knowing whether 
the particles that have arrived may not be the ones we have started with? Theo¬ 
retically, this also translates into a problem for the quantum operators in which 
probability distributions are to be derived from. 

In |ll2j , the one-particle interpretation of relativistic quantum mechanics was 
reinforced by interpreting the wavefunction of a single free particle to be one 
of the two independent (’positive’ or ’negative’) solutions of the Klein-Gordon 
equation (for spin-0 particles) or the Dirac equation (for spin-1/2 particles). 
One-particle analysis would then involve restricting ourselves into a Hilbert 
(sub)space spanned by the positive or negative states and working with true 
one-particle quantum operators which do not map these two subspaces into 
each other. They have calculated several one-particle operators for, say, energy 
(the Hamiltonian), momentum, and position. However, a time operator was not 
constructed implying the parametric role of time in the one-particle formalism. 
This leaves us a certain sense of disparity with regards to the role of time in rela¬ 
tivistic quantum mechanics. That is, in special relativity, time is a quantity that 
is intrinsic to the system in question. Moreover, space and time are intertwined 
so that they are treated on the same footing. However, in the earlier days of 
quantum mechanics, time is not treated as a physical observable quantity but 
only an external and absolute parameter marking the evolution of a system. 
This is in contrast with the observable status of space. The powerful formalism 
of quantum field theory (QFT), which states that particles and antiparticles are 
just excitations of a quantum field subject to fluctuations, demotes the role of 
space so that time and space are on equal footing. They both, however, are just 
parameters labeling the quantum field operators. 

Going in a different direction, there are studies in both non-relativistic 
and relativistic quantum mechanics promoting the role of time to an observ¬ 
able [HE]. Particularly, it was calculated that the TOA operator for a non- 
relativistic particle with mass too is —tooT-iu where the Tm,n’s form 

a complete and linearly independent set called the Bender-Dunne operators 
. In [E] , a proper time operator was found for a relativistic electron and 
was used to study some properties of several position operators. And in [E], 
a method for constructing TOA probabilities valid for any experimental setup 
(including relativistic systems with interactions described by QFT) was devel¬ 
oped. They associated the TOA of a particle as the time instant when there is 
a transition in the degrees of freedom of the detector. While in m, a relativis¬ 
tic particle moving in (5-1-1) dimensions, where two of the spatial dimensions 
are compactified into a torus, was studied. They constructed a discrete physical 
time based on a quasi-local invariant observable. It was then found to be related 
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to the proper time on average. Another study [ms] is on the construction and 
properties of a self-adjoint relativistic TO A operator for a free spin-1/2 parti¬ 
cle. They found that the self-adjointness was due to the existence of particle 
and antiparticle solutions of the Dirac Equation. An important concept that 
has risen from these studies is the idea of supraquantization [7]. It is the idea 
of constructing quantum observables without quantization. This is so that we 
can break from the circularity of promoting a classical quantity into a quantum 
operator then checking if it indeed reduces to the said quantity in the classical 
limit. That is, if classical mechanics can be thought of as only a limiting case of 
quantum mechanics, then quantum mechanics must be autonomous from clas¬ 
sical mechanics. Quantum operators should then be constructed from quantum 
mechanical principles with the added condition that the corresponding classical 
quantities indeed arise in the appropriate limit. TOA operators are then con¬ 
structed as being canonically conjugate with the system Hamiltonian. For a free 
particle, the TOA operator is indeed —moTLi.i [Sill]. For potentials of higher 
order than quadratic in position, the usual quantization of the classical TOA 
does not produce an operator canonically conjugate with the Hamiltonian. Such 
instances are the shortcomings of quantization. However, we can construct a 
TOA operator that is conjugate with the Hamiltonian (from a quantum princi¬ 
ple) and see that it would reduce to the correct classical TOA [5]I5] . Likewise, we 
should then be able to construct a TOA operator for a relativistic particle that 
is conjugate with the system Hamiltonian without referring to any quantization 
rule. As stated earlier, this operator should be made into a true one-particle 
operator since we are in the relativistic quantum regime. 

This study is then a follow-up to [22]. That is, we will construct a relativistic 
one-particle TOA operator for a free spin-1/2 particle in one spatial dimension. 
Here, we will use similar methods done in [mill], but for a free spin-1/2 
particle. In Section [2l we briefly review the one particle formalism of the Dirac 
Equation and define the system Hamiltonian in the process. In Section [S] we 
construct a TOA operator using a method similar to [22] and put the resulting 
operator into a representation where the Hamiltonian is diagonal (the so-called 
<I>-representation). In Section 0] we impose that the TOA operator be invariant 
under parity inversion. And lastly in Section [S] we conclude. 

2 One-Particle Interpretation of the Dirac Equa¬ 
tion 

In (1-1-1) dimensions, the free Dirac Equation may be written as 

ih-— = —ihcA— -h (1) 

at ox 

where A and B are 2x2 hermitian matrices satisfying {A,B} — AB -|- BA = 0 
and A^ = B^ = (q ?), and = {ipi 4 ’ 2 )'^ is a two-component column matrix 
describing a free spin-1/2 relativistic particle with mass mo. The corresponding 
probability density is positive definite and can be readily constructed as p = 
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= (^ip* ip 2 ){i’i ^” 2 )^ = IV'iP + IV' 2 p- We then may define the Hamiltonian 
of the system as 

d 

TL^ = —ihcA ——h moc^B (2) 

ax 

which we call as being in the Schrodinger representation ('^-representation) as 
in mug. Acting twice on 'h, we have 


- A 


92 ^ 

'W 


= ^2 ^ = -h^c^ 


^2^- 

i9a;2 


mnc'^'h 


(3) 


which means that each component of '1' satisfies the Klein-Gordon equation 
separately. 

We also see that the solutions of Eq © are of the form 'h ^ exp{ipx/h — 
iEt/H) which tells us that the eigenvalules of 'Hip are E = XEp = ±Ep = 
±\/p^c2 -|- TOqC^. The Dirac Equation can then be taken as having two types of 
solutions characterized by A, similar with the Klein-Gordon case. The difference 
is that, given a momentum p, the particle indeed can have an energy +Ep or 
-Ep 0 . The interpretation provided by Dirac is that all of the negative energy 
states are occupied whereas the positive energy states are not. The transition of 
a particle from the ’sea of negative states’ to a positive state leaves a hole in this 
sea which we now interpret as an antiparticle. The Dirac sea however necessarily 
requires an infinite number of particles so that a one-particle interpretation 
may not stand a rigorous treatment. Still, we insist on the interpretation by 
considering one-particle operators mm- First, we see that the Hamiltonian 
Hip = cAp + moc^B (in momentum representation) can be diagonalized. Since 
A = and B = PjCTj, where 



we can use 


[/ = = t/-i = 


yj2Ep{Ep + aspc Psmoc^) 


{EpT^ -l- Hip) 


(4) 


to show that 


d 


d 

zfi—$ = H<p$ 
at 


(T3.Ep<^. 


(5) 


We see then explicitly that the eigenvalues of the Hamiltonian are EEp. 

It is interesting to note that the time evolution of states generated by the 
diagonalized Dirac Hamiltonian, Eq ([5]) , is of the same form as the time evolution 
generated by the diagonalized Klein-Gordon Hamiltonian [Tl[2l[^. Now, if we 
consider restricting our wavefunctions to be of the form $+ = (0+ 0)^ or 
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$_ = (0 then <i)± has energy ±Ep. Moreover, Eq ([5]) states that the time 
evolution of $+ is independent of <f>_ and vice-versa. That is, the states do not 
mix as time evolves. We may then consider as a one-particle wavefunction. 
True one-particle operators then should not map a positive state $+ into a 
negative state $_ and vice-versa. So if an operator in the ^'-representation 
is a true one-particle operator, then A<s> = UA^iiU~^ in the so-called Feshbach- 
Villars representation ($-representation) should be diagonal [UIl]- An example 
of a one-particle operator is = cr^Ep. These types of operators do not mix 
the positive and negative states. 

The Hilbert space of states % may then be seen as being split into two 
subspaces H+ and so that % = %+ 0 "H-. To be more specific, since we 
may write the inner product in H as < >= f = / 4'|C7'l'C/4'2dp = 

f ^\^ 2 dp =< $i |$2 >, then 


Ux 



( 0(A) 
1^0(-A) 


'('A 


< $a|$a >= 



\(j)x\^dp < oo 


with A = ±1. If we, say, take $+ = (1 as describing a particle, then 

analysis of that same particle means that we are restricting ourselves in the 
state subspace T-L^ and are working with operators which do not mix the states 
of 'H+ and %- 


3 Constructing 7^ 

In this section we attempt to construct a Time of Arrival operator canonically 

conjugate with the system Hamiltonian Hip so that 


Hip, 


'k = ill'll 


( 6 ) 


for any two-element column matrix 'k, and the Hamiltonian is given hy T-L^, = 
cAp + moAB. We proceed by expressing 7i in terms of the complete linearly 
independent set of basis operators (so-called Bender-Dunne operators) 

[201I1I], 


1 


n 


n\ 




^kpm^n-k 


(7) 


Explicitly, 


^ ( Crn^n^m^n 


( 8 ) 


where the Cm,n^ just constant 2x2 matrices. Also, since the CTj’s form 
a complete linearly independent set of 2 x 2 matrices, any 2x2 matrix can 
then be written as a linear sum of the Uj’s. Specifically, we can write Cm,n = 
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S?=o where 7 ™’" are just unknown scalars. Substituting Eq ([5]) into Eq 

([ 6 |, we have 




^ ^ .Ap, Crn,n^m,n “t TTloC B^ Cm^nTm^r. 
m,n 

^ ^ (^^C7n,nP'^m,n ^m,n^^m,nP^ H“ ^77-0^ [-S, Cm,n\ ^m,n^ ^ 


m,n ^ 

= Cm-l,n] — 2 


ihcn , 


m+l,n 2 {^5 C^m,n}7j.n,n— 1 “t“ TTIqC [B^ C^m,n] ^m,r, 


ihc{n + 1 ) r . ^ , 2 


— ^ ^ ^h(To(jm,0^n,0^m,nlk 


C^m,n+l} “t“ TUqC [B^ C^m,n] j 

(9) 


where in the fourth line, we shifted indices for each term and in the last line, we 
assert the equality of the right hand side of Eq ([5]). Then, the coefficients of the 
Tm.n’s should be equal. Additionally, since the aj’s are linearly independent, 
their coefficients should also be equal. These would restrict the values of the 
coefficients 7 ™’"’s. 

Strictly speaking, the coefficients 7 ™’*^ would remain arbitrary. However, we 
wish to take as many coefficients as possible to vanish. That is, we wish to take 
the minimal solution so that we can set 7 ™’° = 0. It can be shown that the 
only non-vanishing coefficients are 23 }“ 2 3 } ~ 

Explicitly, we then have 

Tq/ = C-iqTLin + Co^iTq^i 

= -mo(/ 3 icri -I- P 2 cr 2 + P 30 ' 3 )T-i,i - -(aiCTi + 020-2 + aao-ajToq 

c 

= —moBT-ii - ATo,i 

( 10 ) 


so that in momentum representation (specifically, 'k — p representation), its 
action on a column matrix d/ (p) is 

( 7 *) (p) = -mosf (i|»M + I Qt(p))) - ( 11 ) 

Much like in the spin-0 case, we transform our operators to the $ — p rep¬ 
resentation. This transformation is made through 7$$ = for some 

arbitrary column matrix $. The unitary transformation [/ diagonalizes the 
Hamiltonian so that 77$$ = U'Hq,U~^^ = E^az^- Its eigenstates are then 
explicitly separated into positive and negative states. This makes one-particle 
analysis more straightforward. Explicitly, U is given by Eq (|4]) so that we can 
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calculate 


= -moUBT.i - -UATq lU 

’ c 


-1 


$ 


= -moUBU-^T_i 1$ - -UAU-^To - ihU 

’ c 


A moB\ dU ^ ^ 


P 


= -moUBU 


-1 


ihd^ ih 


ih - 19$ 

= — ^un^u-^ 


p dp 

ihU 


- 1 - -uAu-\h^ - 

C op pc^ op 


pc^ 

- F 

— 9 -^P^3 

pc^ 


where we can let 

- ^ ih ^ 

T$ =- 

pc^ dp 


ih 


dp 

9$ mg'’^ 


V ■ 


1 


9p pc 
dU-^ 




pc^ dp 


., rUnC 

+ P's*!* - 


(ai/32 - Q!2/3 i) 


$ 

hmoc 


dp 2p'^Ep 


2 2 

■fc ^ 

■ ih ^ g3$ 


(i?p + peas + moc^Ps) 2p 


$ 


( 12 ) 


2p2£;p 


= -:T7r 2p + 2 


2R 




9$ TOnC^ 
- 


dp 


2Ep 

1 


9<I> ih 


dp 


= -TTTTO'sT"!,!^ - 


2R 


= -Air^3 ihp^ + + 2tooC )a3 


2Eri 


2E, 

+ § I ^3r-,,.<5 


19$ ift 


p dp 2p‘^ 


ro$ = - 


(ai/32 — Q!2/3 i) hmoc 

{Ep + peas + moc^Ps) 2p 


$ 


(ai/32 - a2/3i) 


hmr 


{Ep/c^ + pas/c +molds) 2pc 


$ 


(13) 


We see that T has the same action for the spin-0 case. It was shown that in the 
non-relativistic limit, T reduces to —moT-i^i which is just the TOA operator 
for the free non-relativistic particle. Also in the non-relativistic limit, it can be 
shown that To vanishes. 


4 Imposing Parity Inversion Symmetry 

Consider the one-particle parity operator 11, whose action on a column matrix 
$(p) in $ — p representation is n$(p) = $(—p). This allows us to calculate the 
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following 


nr$(p,t) 

tn$(p,f) 


ih 

- yEpas 


ih 

- nEpas 

pc^ 


dp 

dH-p,t) 

dp 


+ ih 
+ ih 


TOgC^ 

2p^Ep 

m-gC^ 

2p^Ep 


a3<^{-p,t) 


a3<^{-p,t) 


n,r 


$(p,f) = 0 


( 14 ) 


where, df{x)/dx\x=-p = —df{—p)/dp. This implies that T is invariant under 
parity inversions, i.e. 117^11“^ = T. Similarly, 


IlfQ^{p,t) 

foll^{p,t) 





(ai/^2 - a2l3i) 

{Ep - peas + mgc^/^a) 2p 

(ai/32 - Q2/3 i) .n 

{Ep + pca3 + TOoc^/Ss) 2p 
2(^p + wgc^/33)(ai/32 - a2/3i) ^^gc 

+ TOgc 2 / 33)2 - (pcaa)^ 2 p 


(15) 


We see then that To is not invariant under parity inversions. 

If we impose this symmetry on our TOA operator for free spin-1/2 particles, 
then To must be thrown away. Since it can be shown that [H^, To] = 0, we can 
interpret To as a term not contributing to the conjugacy of the TOA operator 
with the Hamiltonian Actually, we can add any operator that commutes 
with to one which is conjugate to 7/$, suggesting that there are many 
solutions to Eq (jS]). These operators then represent the many characteristics, 
roles, and physical interpretations of a time observable [Mill]- We are certain 
that the remaining term T” represents a TOA observable since its action is the 
same with the TOA operator from and since the time evolution of states, Eq 
®, is similar to that of the Klein-Gordon case [ijiuiia, the dynamic behavior 
of the eigenfunctions would also be similar to the Klein-Gordon case [35]. 

As an example, the time evolution of the probability density p = = 4)1$ 

(in configuration space) of a typical eigenfunction corresponding to the 

eigenvalue r = 1.0 is given in the contour plots [33]. It can be seen that in 
the contour plot, the probability density peaks about the origin at the 

time of the eigenvalue. Whereas in the P^~'> contour plot, the two peaks are 
closest to the origin at the time of the eigenvalue, but with the probability 
always vanishing at the origin. We call the states with similar behaviors as 
demonstrated in p(+) as a non-nodal eigenfunction while the states with 

similar behaviors as in as a nodal eigenfunction The localization of 

the states at the origin at their corresponding eigenvalues is then interpreted 
as the particle arriving at the origin at the time of the eigenvalue. That is, the 














Position X Position x 


Figure 1: Time Evolution of the Probability Density with t = 1.0 


eigenfunctions represent states of definite arrival time at the origin. The specific 
difference would be the inner products defined on the vector spaces on which 
they act. Hence, a one-particle TOA operator for a free spin-1/2 relativistic 
particle is given by T whose action is the same as the TOA operator for the 
spin-0 case. 


5 Conclusion 

We have constructed the minimal operator solution Tip that is conjugate with the 
Hamiltonian for a relativistic spin-1/2 particle. It is a true one-particle operator 
since it does not map positive states into negative states and vice-versa. That is, 
it is diagonal in the <I>-representation. We see that the action of the said operator 
7i = UThU~^ =T + To is the same as the one derived from [22] apart from an 
additional operator term To which commutes with the Hamiltonian. But upon 
imposing that a TOA operator should be invariant upon parity inversion, we 
see that To breaks this symmetry. We can throw away this term so that our 
relativistic one-particle TOA operator for a free spin-1/2 particle is T. (Any 
other implications of To, if any, may be studied later.) Since the calculated 
actions of our TOA operators for the spin-0 and spin-1/2 free particles are the 
same, and the time evolution of states are the same, we expect similar properties. 
As an example, they should both reduce to the non-relativistic TOA operator 
—tuqT-i i, and the form and characteristics of their eigenfunctions must be the 
same. The difference would be in the inner products defined on the Hilbert 
spaces they are acting upon. Then in this relativistic quantum regime, time is 
an observable represented by an operator. More importantly, it was constructed 
by imposing a quantum principle and hence, without resorting to quantization. 

One of the motivations of studying time in relativistic quantum mechanics 
is to lessen the inconsistencies between the notions of time in special relativity 
and quantum mechanics. The hope is that, from these studies, we may gain 
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further insight that can help us mend the marriage of general relativity and 
quantum mechanics. The reconciliation of which is very much harder. Perhaps 
that insight is supraquantization. Consider the standard canonical quantum 
general relativity approach for example, spacetime is split into space and time 
so that the classical equations describing the system can be cast into Hamilto¬ 
nian form |231I24] . This Hamiltonian however, is actually constrained to vanish 
via the equations of motion. The quantum version of this is supposedly to im¬ 
pose the vanishing of the quantized Hamiltonian as a condition on the physical 
states. This results into the Wheeler-DeWitt equation - similar in form to the 
Schrodinger equation but with the time derivative vanishing or ’missing’. This 
vanishing leads to one facet of the problem of time in quantum gravity [25] (Is 
the universe ’frozen’?. Is time an illusion?). Perhaps this is where supraquan¬ 
tization can come into play. Perhaps there is a more ’fundamental’ way to 
construct operators, like the Hamiltonian stated earlier, by using some rele¬ 
vant quantum first principle(s) so that we won’t need quantization. Of course, 
even determining the appropriate first principles would require further studies 
to plant these speculations on firmer and more rigorous grounds. Nevertheless, 
these studies may lead us closer on formulating a theory of quantum spacetime. 


A Calculation Details of the Construction of 7^ 


We consider the commutator between the position and momentum operators q 
and p, respectively, [q,p] = ih. We calculate the following: 


= -ihq^-^ + {q’^-^ [p, q] + [p, q’^-^] q)q 

= -2ihq^~^ + f 


ihkf-^ 


(16) 


Now, consider the complete and linearly independent set of Bender-Dunne op¬ 
erators denoted by Tm,n- 




( 17 ) 
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and calculate the following: 


oo 

-Y 

2” ^ 

n! 

k\(n — 

k)\ 

ip<tW 

Ti ^n — k 

1 

2n ^ 

n! 

k\{n — 

k)\ 

■Yp- 

ihkq'^-^)p'^q'^-^ 


ihn 




^m+l,n 

~ ~ 

Tm 

,n—l 


1 

2n ^ 

n! 

k\(n — 

ky. 


- oo 

2” ^ 

n! 

k\(n — 

k)\ 

+ ihkf-^) 


ihn 




^m+l,n 


Tm 

,n—l 



(18) 


We also consider a pair of 2 x 2 matrices A = ^ ~ ^J=o 

satisfying the anti-commutation relation {A, B} = 0 and A^ = = ao where 

the set of aj ’s form a complete linearly independent set of 2 x 2 matrices, denoted 

by 



Note that ao is just the identity matrix. Some properties of these matrices 
include axay = iaz where, {x,y,z) are just cyclic permutations of (1,2,3), and 
ajaj = (To for j = 0,1, 2, 3. The condition that A and B anti-commute and that 
both of their squares equal to the identity matrix imply that we should have the 
following conditions «□ = /^o = 0, oijPj = 0, and P] = ^- 

We also restrict the a^’s and /3j’s to be real so that A and B are self-adjoint. The 
Hamiltonian = cAp + moc^B is then also self-adjoint. Using the properties 
of the (Tj ’s, we can calculate some relevant commutation and anti-commutation 
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relations of A and B with Cm,n = Sj=o Specifically, 

[A, ai] + 72"*-'’”[71, a2] + ir^'^[A, ag] 

_ 2ia2cr3 + 2103 ( 72 ) + 7™~^’”(2iai(j3 — 2m3(7i) 

+ 7r~^’"(~2ioi(72 + 2ia2(Ti) 

= (^(73 “2 - 72 «3)cti + (7i as-ls “ 1)^2 

+ (72 ai - 7i 02 ) 0-3 j 

[B, Cm,nl = 2Z ((73"*’”^2 - 72”^’”/33)ai + (ir^S - iT'"^l)o2 
+ (72™’”/3i-7r’"/32)a3) 

{71, = 7o"*’”+'{^, txo} + 7r’”+'{A, (Til + 72™’"+'!^, ^ 2 } + 73™’"^'{^, ^^ 3 } 

_ 2-)/™’”+i(aicri + 020-2 + 030 - 3 ) + 27™’”'''^ai(To 

I n m,n+l , n m,n+l 

+ 272 020-0 + 273 O 3 O -0 

n/ m,n+l . m,n+l . m,n+l \ 

= 2(73 Oi + 72 ’ 02 + 73 ’ 03 )(To 

I n m,n+l . n m,n+l . n m,n+l 

+ 27 o oi(Ti + 27 o 020-2 + 27 ^ 03(73 

(19) 


From Eq Q, we then arrive to four equations restricting the values of the 


m.rif 

7,- s. 


-'-fin + 1)2 (7r’"+'oi + 72 ™’"+'o 2 + 73"*’"+'a3) = ihSm.oSu.o (20) 

n-/ m —l,n m —l,n \ i^C , i\o m,n+l 

c2z(73 ’ 02 - 72 “ 3 ) - — (^^ + 1)27o “i 

+moc22*(73”^’"^2 - 72”^’”/33) = 0 (21) 

o-/ m —l,n m —l,n \ ihc , > m,n+l 

c2i(73 ’03-73 “1) - —+ 1)27o «2 

+moc22z(7r^’”^3 - = 0 (22) 

o ■ / m —l,n m —l,n \ i^C , i\o m,n+l 

c2i(72 ’01-73 “2) - — (?T- + 1)27o 03 

+moc22z(72”^’”^i - 7r’”/?2) = 0 (23) 


Using Eqs (I21II23|) . to cancel their respective third terms, we arrive at 

7r“'’"(/3203 - / 3302 ) + 72"*-'’"(/330i - / 3 i 03 ) + 73™“'’”(/3l«2 - /320 i) 
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Using again Eqs (i^rra . this time to cancel their respective first terms, we 
arrive at 


-moc(7™’"(/32Q;3 - ^ 0 . 2 ) + 7™’"'(/33ai - , 8103 ) + 7r’"(/3i“2 - hai)) 

h 

- -(n + 0 

(25) 

so that Eq (l24l) and Eq (1251) imply that = 0. Eor n ^ 0, Eqs (I20II23I1 can 

then be written as 

/ m,n I m.n . m.n \ c r /r\n\ 

-cn ( 7 i Oi + 72 ’ 02+73 “a) = Om,oOn,l (26) 

m — l,n I m.n o m—l,n . m.n n /n^\ 

73 ’ 02 + moC 73 ’ 82 = 72 03 + TO 0 C 72 83 (27) 

m —l.n I m.n 0 m—l,n , m.n n /no\ 

73 ’ 03 + moC 7 i ’ = ’ oi + TO 0 C 73 81 (28) 

m — l.n I m.n 0 m—l,n . m.n n /nn\ 

72 ’ oi + moC72 81 = 7i 02 + tooC7i 82 ( 29) 


Appropriately combining Eqs (|?7ll and ([28]), (l27|l and (l29|l . and (l2^ and (l29ll . 
respectively, we arrive at 


/ m—l.n I m—l.n \ . / m.n . m.n \n 

03(71 ai+72 a2 )+TOoc( 7i oi + 72 ’ 02)83 


l-o| 


m—l.n/ 2 I - 2 ^v^i m.n/ , 

= 73 i£ 4 +^ 2 ) + W0C73 jai8i-t^02/: 

/ m—l.n I m—l.n \ . / m.n . m.n \ 

02(71 ai+73 03 )+TOoc(7i 01+73 ’ “3)82 


1 - o| 


7^ ^’"i07+“+f|r+’7roC7r’"j3bS^^ 


/ m—l.n I m—l.n \ . / m.n . m.n \ n 

ai(72 “2+73 03)+TOoC(72 02 + 73 ’ “3)81 

1 — of 


m—l.n/ 2 I I m.n/ 

= 7i ia|+^) + woC7i (c^ 


—0383 


—0282 


—0181 


and using Eq (gnj), we get n 7 ™:|^^ 2 , 3 } = -Sm,o^n,iaj/c- moSm,-iSn,il3j. Note 
however, that there are some arbitrary coefficients 7 j(L’|o 123 } which do not 
necessarily vanish. We wish to take the minimal solution so that we take to zero 
as many 7 ^’^’s as possible. That is, we set the arbitrary coefficients 7 ™’° = 0. 
The only non-vanishing coefhcients are then 23 }“ ^ ^J={^ 23 }“ 

-TOo8i- 
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B Calculation Details of the Transformation 7$ = 

u%u-^ 


We first take the derivative of Eq (jd]) 
dU~^ _ (Ep + oapc +( pc^ 


dp 


V2 


. 2E. 


, 3/2 


cA pc^ ~ 

^3 d- 77^ -EVTT^'I' 


2Ep 


{Ep + aapc + Psmoc^) (pc^ 


Er, 


- a^c 


, 5/2 
p 

EpCr^ + 


E. 


,1/2 


2 V 2 

_ (Ep + aspc + ' 

2Epy^2E^ 

+ (pc^as + 2cEpA - ) (Ep + a^pc + dsmoc^) 


-{pc^ + cEpaz){Ep(Jz + "H®) 


_ (Ep + Q!3pc + Pampc^) 

‘2Ep \J 2Ep 


[{p^c^a^ - cEpUs + pmoc'^f33)a3 


+2{cEp +pc^Epa3 + moc^Epj33)A - (^pc^ + ^ + cEpOa + 


Ep 

(30) 


Remembering that %% = E^a^ = {p^<^ + TO§c^)croi we calculate 


7^4, i9C/ ^ _ (Ep + oapc + j53m3c^) 2 
pc^ dp 2Ep^2Ep 

+ 2 (S + e ,„3 + !!!!!£aA' 

\ pc P 


la 'Cn\(?OL3 


'moc^d3 - 


'H^O'3 


{pcao + ttloC^BA) 


- f 2Ep +pcEpa3 H-+ moc^Ep^a j ao 


{Ep + a3pc + j33moc^) 2 


+2 


2Ep\/2Ep 
/ mocE,2 


2 a "loC^«3 
toqc da- 


{pcAu3 + TOoc^Eaa) 


'p TOoC^^p/^3 \ / 2 ^ mnC^oa^ 

— + moC^Epa3 H-—^— j EA + ( moC^/3a-^-j Epcro 


( 31 ) 
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Afterwards, we also calculate 


moBU 


so that 
moBU 


{Ep + aspc + iSsmoc^) ""^nioEpB ^ ~ 

— = - — — 7^ -5 —[Ep + aspc + Psuioc ){Epaz + %^) 

^-t^p yj P 

^/moE^ ^ mocEpas ^ mlc^Epfis'^ 


_ {Ep + aspc + 2 


p^ 


P 


‘^Ep^j2Ep 

m^El mocEpas mlc^Epl3:i\ ^ ^ 2 ^ 

- 5 -^ H- - -^- o — {pcBA + TOoc (To) 

'n^ n I 


pz 


Bas 


p 

\ 

{Ep + aspc + l3^moc^)~i 

‘2Ep^j2Ep 


pZ 

/ 

/ mocEl 




P 


- moC^Epas 


mlc^EpPs 


BA 


’m-pC^E^ mlc^Epa^ mlc'^EpP^ 


CTo 


\ " " ' / 

/ moE{^ rrinC^as 2 4o ^ \ „ 

-5^ + pnioc^as H-^ + TOoC^/^s H-5-/33 -8(73 

\ P P P J 


(32) 


we can subtract them 


1 n^, dU ^ {Ep + aspc + l33nioc^) = r/ 2 4 3 o ^3 

- -- [<"■»' “» - 


+ 


/ mo-Ep / 3 „TOoC®\ WqC® \ 

( ^ + ( pmoc® + 2——\ a 3 + 1 - 8(73 

^ ^mlc^Epa^i ^ ( mj^c^ ^ 2 ^ i:. « ^ _ 

1 ^2 ^ ^ + (^-^5 mpc j EpPs j (70 


+ 


I mpcEp , _ 27^ , '" -o^ ^ph'i I 




P 


2p , mlc^EpPs 
- nioc EpUs 


P 


(33) 
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Lastly, we muliply the result by U. 


/moBU-^ 

V 2p2 


'Hqf dU {Epas + pcA + moc^B) ^ moBU ^ T-L^, dU ^ 

pc^ dp ) y/2Ep{Ep + a^pc + jS^m^c'^) ^ 2p2 pc^ dp 


(Ep + aspc + Psmoc^y 

4E2 


[(pniQC^as - p^moc^pyA'^as 


5^8 


( TUqC E niQC^ 2 2,0 2 4 ^ I O I d 2 

H- (p C +2mgC jas H- 5 —ps B as 


\ P'^ P 

/mgC^Llp Ic^EpUs mo 


moC 




+ 2 


( mlc^El 


. / 2 4 2 2\ 7:^2 /O 

H- Y -pc) fjaao 


+ 2mgC^Llpa3 + 


- pmgc^ EpPs Aao 


+ tjiqC^ EpOts “t“ 


3^5 p o \ 
]L -g/pP3 I 


BAB 


p - ■ p J 

+ {mlc^Epas - pm^c^EpPs)asAas 


ml^El ^^mle'EpUs , (mlc^ _2„4 


pz 


P 


pZ 


- m^c EpPs Berg 


+ [moyEp + pmoc^Epas + mlc^EpPs) A^B 


moEl ( mle'Er, 

prriQC Ep + 2- 


\ P' 

( mocE^ 


P 


as + 


niQC^Ep 


P^ 


Ps asBas 


V P 

( mscEl 


„ „ „ wlc^E'iPs\ 

^ + moc^Epas H-^— j ugAB 


P 


4 7 
m^c 


—- + (p'^nioC + 2tooC®) as H- -—Ps ABas 

P P J 

3„6„, „™2„5/ 


+ {rriQC^as - pm^CPs)BAas 

{Ep + aspc + Psmoc^y"^ 


4B2 


moE^ 

p2 


pmoc Ep + 2- 


[{mlc'^^Epas - pmoc^EpPs){A + crgAcrg) 
m^e’E, 


tQ'-' J^p 

P 




^ mgc^B^aa , ^mg ^ 2 4 2 2 ^ r .20 \ 

2- 5 —^ + 4--h 2— (TOgC -pc) E Ps as 

pZ p pZ \ / IJ j 


/ mocE^ 

^ \ P 


„ „ wlc^ElPs\ , , , , 

+ moc^Epas H-^— ) (agAB + ABas) 


P 


(34) 
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{Ep + aspc + /?3moc^) 
4^2 


2\-2 


mlc^El ^mlc^E'^as rr^cpi]p3s 




P 


+2myE,al + + 4l^fc!:5«A'| 


P 


0-3 


/ niocE^ 

\ P 


+ nioEEpas + 




,2^3 ^2/ 

~P ) 


{Ep + aspc + Psmoc^y 

4^2 


2EpmQC^ 

p2 


{Ep + aspc + Psmoc ) as 


rriQcEp „ 

-{Ep + aspc + / 33 TO 0 C )2i(ai/32 - a2/3i)o-o 


P 


2^2 


rriQC 

2p‘^Ej 


-as 


moc 


2p{Ep + a3pc + Psmocy 


i{aiP2 - a2l3i)ao 


(35) 


Note that {^,-B} = 0 and = B'^ = ag, where A = X)j=i ^ = 
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